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At very high densities, electrons react with protons to form neutron rich matter. This 
04 ' material is central to many fundamental questions in nuclear physics and astrophysics. More- 

over, neutron rich matter is being studied with an extraordinary variety of new tools such 
as the Facility for Rare Isotope Beams (FRIB) and the Laser Interferometer Gravitational 
Wave Observatory (LIGO). We describe the Lead Radius Experiment (PREX) that uses 
, \ parity violating electron scattering to measure the neutron radius in 208 Pb. This has im- 

^ . portant implications for neutron stars and their crusts. We discuss X-ray observations of 

neutron star radii. These also have important implications for neutron rich matter. Gravi- 
tational waves (GW) open a new window on neutron rich matter. They come from sources 
such as neutron star mergers, rotating neutron star mountains, and collective r-mode oscil- 
(N. lations. Using large scale molecular dynamics simulations, we find neutron star crust to be 

very strong. It can support mountains on rotating neutron stars large enough to generate 
detectable gravitational waves. We believe that combing astronomical observations using 
photons, GW, and neutrinos, with laboratory experiments on nuclei, heavy ion collisions, 
and radioactive beams will fundamentally advance our knowledge of compact objects in the 
' heavens, the dense phases of QCD, the origin of the elements, and of neutron rich matter. 

§1. Introduction 

Multi-messenger astronomy observes matter under extreme conditions. In this 
paper we describe how electromagnetic, gravitational wave, and neutrino astronomy, 
along with laboratory experiments, provide complimentary information on neutron 
rich matter. Compress almost anything to very high densities and electrons react 
with protons to form neutron rich matter. This material is at the heart of many 
fundamental questions in Nuclear Physics and Astrophysics. 
| • What are the high density phases of QCD? 

• Where did the chemical elements come from? 

• What is the structure of many compact and energetic objects in the heavens, 
and what determines their electromagnetic, neutrino, and gravitational-wave 
radiations? 

Furthermore, neutron rich matter is being studied with an extraordinary variety 
of new tools such as the Facility for Rare Isotope Beams (FRIB), a heavy ion ac- 
celerator to be built at Michigan State University,^ and the Laser Interferometer 
Gravitational Wave Observatory (LIGO)P^ Indeed there are many, qualitatively dif- 
ferent, probes of neutron rich matter including precision laboratory measurements 
on stable nuclei and experiments with neutron rich radioactive beams. While as- 
trophysical observations probe neutron rich matter with electromagnetic radiation, 
neutrinos, and gravitational waves. In this paper we give brief examples of how 
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neutron rich matter is being studied with these extraordinarily different probes. 

We are interested in neutron rich matter over a tremendous range of densities 
and temperatures were it can be a gas, a liquid, a solid, a plasma, a liquid crystal, 
a superconductor, a superfluid, a color superconductor, etc. Neutron rich matter 
is a remarkably versatile material. The liquid crystal phases are known as nuclear 
pasta and arise because of coulomb frustrationP'® Pasta is expected at the base 
of the crust in a neutron star and can involve complex shapes such as long rods 
( "spaghetti" ) or flat plates ( "lasagna" ) . Neutrinos in core collapse supernovae may 
scatter coherently from these shapes (neutrino pasta scattering) because the shapes 
have sizes comparable to the neutrino wavelength.^ 

In this paper we focus on some of the simpler gas, solid, and liquid phases of 
neutron rich matter. In Section [2] we describe a precision laboratory experiment 
called PREX to measure the neutron radius of 208 Pb. Nuclei are liquid drops, so 
PREX and many other laboratory experiments probe the liquid phase of neutron rich 
matter. In astrophysics, electromagnetic, gravitational wave, and neutrino probes 
can observe different phases of neutron rich matter because the probes have very 
different mean free paths. In Section [3] we describe electromagnetic observations 
of neutron star radii. In Section [5] we discuss gravitational waves from neutron star 
mergers that are produced by the energetic motions of dense liquid phase neutron rich 
matter. In addition, continuous gravitational waves can be produced by "mountains" 
of solid neutron rich matter on rapidly rotating stars. We conclude in Section [5j 

§2. Laboratory probes of neutron rich matter 

Neutron rich matter can be studied in the laboratory. Hot and or dense matter 
can be formed in heavy ion collisions, while more neutron rich conditions can be ac- 
cessed with radioactive beams. In addition precise experiments are possible on stable 
neutron rich nuclei. We give one example, the Lead Radius Experiment (PREX) 61 
accurately measures the neutron radius in 208 Pb with parity violating electron scat- 
tering.^ This has many implications for nuclear structure, astrophysics, atomic 
parity violation, and low energy tests of the standard model. 

2.1. Introduction to neutron densities and neutron radii 

Nuclear charge densities have been accurately measured with electron scattering 
and have become our picture of the atomic nucleus, see for example refP These 
measurements have had an enormous impact. In contrast, our knowledge of neutron 
densities comes primarily from hadron scattering experiments involving for exam- 
ple pions,® protons E3>[ID>[I3 Q r antiprotonsP^'^ See also refP^ for a beautiful 
measurement of the dipole polarizability. However, the interpretation of hadron 
scattering experiments is model dependent because of uncertainties in the strong 
interactions. 

Parity violating electron scattering provides a model independent probe of neu- 
tron densities that is free from most strong interaction uncertainties. This is because 
the weak charge of a neutron is much larger than that of a protonPSJ Therefore the 
Z° boson, that carries the weak force, couples primarily to neutrons. In Born ap- 
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proximation, the parity violating asymmetry A pv , the fractional difference in cross 
sections for positive and negative helicity electrons, is proportional to the weak form 
factor. This is very close to the Fourier transform of the neutron density. Therefore 
the neutron density can be extracted from an electro-weak measurement.^ 1 Many 
details of a practical parity violating experiment to measure neutron densities have 
been discussed in a long paperP 

The neutron radius of 208 Pb, R n , has important implications for astrophysics. 
There is a strong correlation between R n and the pressure of neutron matter P 
at densities near 0.1 fm -3 (about 2/3 of nuclear density) P^-* 1 A larger P will push 
neutrons out against surface tension and increase R n . Therefore measuring R n con- 
strains the equation of state (EOS) — pressure as a function of density — of neutron 
matter. 

Recently Hebeler et alP^ used chiral perturbation theory to calculate the EOS 
of neutron matter including important contributions from very interesting three neu- 
tron forces. From their EOS, they predict R n — R p = 0.17 ± 0.03 fm. Here R p is the 
known proton radius of 208 Pb. Monte Carlo calculations by Carlson et al. also find 
sensitivity to three neutron forces.^ Therefore, measuring R n provides an impor- 
tant check of fundamental neutron matter calculations, and constrains three neutron 
forces. 

The correlation between R n and the radius of a neutron star r^s is also very 
interesting!-^ In general, a larger R n implies a stiffer EOS, with a larger pressure, 
that will also suggest r^vs is larger. Note that this correlation is between objects 
that differ in size by 18 orders of magnitude from R n ~ 5.5 fm to r^s ~ 10 km. We 
discuss observations of rj^s in Section [3l 

The EOS of neutron matter is closely related to the symmetry energy S. This 
describes how the energy of nuclear matter rises as one goes away from equal num- 
bers of neutrons and protons. There is a strong correlation between R n and the 
density dependence of the symmetry energy dS/dn, with n the baryon density. The 
symmetry energy can be probed in heavy ion collisions.^ For example, dS/dn has 
been extracted from isospin diffusion data 22 ' using a transport model. 

The symmetry energy S helps determine the composition of a neutron star. 
A large S, at high density, implies a large proton fraction Y p that will allow the 
direct URCA process of rapid neutrino cooling. If R n — R p is large, it is likely that 
massive neutron stars will cool quickly by direct URCAPB In addition, the transition 
density from solid neutron star crust to the liquid interior is strongly correlated with 
R -R S3 

lL n ± L p • 

Finally, atomic parity violation (APV) is sensitive to it^pSpSP Parity viola- 
tion involves the overlap of atomic electrons with the weak charge of the nucleus, 
and this is primarily carried by the neutrons. Furthermore, because of relativistic 
effects the electronic wave function can vary rapidly over the nucleus. Therefore, the 
APV signal depends on where the neutrons are and on R n . A future low energy test 
of the standard model may involve the combination of a precise APV experiment 
along with PV electron scattering to constrain R n . Alternatively, measuring APV 
for a range of isotopes can provide information on neutron densities!^ 
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2.2. The Lead Radius Experiment (PREX) 

We now discuss a direct measurement of R n . Parity violation provides a model 
independent probe of neutrons, because the Z° boson couples to the weak charge, 
and the weak charge of a proton = 1 — 4 sin 2 Oy/ ~ 0.05 is much smaller than 
the weak charge of a neutron Q"^ = — 1. Here &w is the weak mixing angle. 

The Lead Radius Experiment (PREX) at Jefferson Laboratory^ measures the 
parity violating asymmetry A pv for elastic electron scattering from 208 Pb. The asym- 
metry A pv is the fractional cross section difference for scattering positive (+), or 
negative (-), helicity electrons, 

dcr I da I 

A dfi I + dQ I ~ /n i\ 

PV ~ do_\ _l dg_\ ■ y Z ' L ) 
dn\+ dtt\~ 

In Born approximation, A pv arrises from the interference of a weak amplitude of 
order the Fermi constant Gf, and an electromagnetic amplitude of order the fine 
structure constant a over the square of the momentum transfer g 2 ™ 

Here the weak form factor Fw (q 2 ) is the Fourier transform of the weak charge density 
pw(r), that is essentially the neutron density, F\y(q 2 ) = J d 3 r sin ^T' pw (r) . Likewise, 
the electromagnetic form factor F c h{q 2 ) is the Fourier transform of the (electromag- 
netic) charge density p c h (r) . This is known from elastic electron scattering.® There- 
fore, measuring Ap V as a function of q allows one to map out the neutron density 
p n ( r )- Note that, for a heavy nucleus, there are important corrections to Eq. 12 21 
from Coulomb distortions. However, these have been calculated exactly by solving 
the Dirac equation for an electron moving in both a Coulomb potential of order 25 
MeV and a weak axial vector potential of order electron voltsPS} Therefore, even 
with Coulomb distortions, one can accurately determine neutron densities. Note that 
this purely electroweak reaction is free from most strong interaction uncertainties. 

The PREX experiment measures A pv for 1.05 GeV electrons elastically scattered 
from 208 Pb at laboratory angles near five degrees. The first measurement yielded 
A pv = 0.656 ± 0.060 (statitistical) ±0.014 (systematic) ppmP From this the rms 
neutron radius R n minus proton radius R p for 208 Pb was found to be R n — R p = 
0.33!q]J fm. See also ref. 29 ' for more details of this analyis. A second PREX 
run is now approved to accumulate more statistics and reach the original goal of 
determining R n to 1% (±0.05 fm). 

In addition to PREX, many other parity violating measurements of neutron 
densities are possible, see for example.^ Measuring R n in 8 Ca is particularly at- 
tractive. First, 48 Ca has a higher experimental figure of merit than 208 Pb. Therefore 
a 48 Ca measurement may take less beam time than for 208 Pb. Not only does Ca 
have a large neutron excess, it is also relatively light. With only 48 nucleons, mi- 
croscopic coupled cluster calculations,^ or no core shell model calculations p^ 1 may 
be feasible for 8 Ca that are presently not feasible for 208 Pb. Note that these mi- 
croscopic calculations may have important contributions from three nucleon forces. 
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This will allow one to make microscopic predictions for the neutron density and re- 
late a measured R n to three nucleon forces and in particular to very interesting three 
neutron forces. 

§3. Electromagnetic observations of neutron star radii 

The structure of a neutron star can be calculated with the Tolman-Oppenheimer- 
Volkoff Equations of General Relativity 33 ' and is completely determined by the equa- 
tion of state of neutron rich matter. The radius of a neutron star depends on the 
pressure of neutron matter at normal nuclear density and above, because the central 
density of a neutron star can be a few or more times that of normal nuclear density. 
A higher pressure will lead to a larger radius. It is important to have both low 
density information on the equation of state from PREX, and high density infor- 
mation from measurements of neutron star radii. This can constrain any possible 
density dependence of the equation of state from an interesting phase transition to 
a possible high density exotic phase such as quark matter, strange matter, or a color 
superconductor. For example, if the 208 Pb radius is relatively large, this shows the 
EOS is stiff at low density (has a high pressure). If at the same time, neutron stars 
have relatively small radii, than the high density EOS is soft with a low pressure. 
This softening of the EOS with density could strongly suggest a phase transition to 
a soft high density exotic phase. 

The radius of a neutron star rjvs can be deduced from X-ray measurements of 
luminosity L and surface temperature T, L = Airr 2 NS asBT ■, with o~sb the Stefan 
Boltzmann constant. Recently Steiner, Lattimer, and Brown have deduced masses 
and radii 3 -^* from combined observations of six neutron stars in two classes: 1) X-ray 
bursts, and 2) neutron stars in globular clusters. They conclude that observations 
favor a stiff high density equation of state that can support neutron stars with a 
maximum mass near 2 M & and that the equation of state is soft at low densities so 
that a 1.4 M & neutron star has a radius near 12 km. They go on to predict that the 
neutron minus proton root mean square radius in 208 Pb will be R n — R p = 0.15 ±0.02 
fm. Note that this is a prediction for a nucleus based on an equation of state 
deduced from X-ray observations of neutron stars. The Steiner et al. paper^ is 
potentially controversial because their results depend on, among other things, the 
model assumed for X-ray bursts. 

§4. Gravitational Waves 

We turn now to gravitational wave observations of neutron rich matter. Albert 
Einstein, almost 100 years ago, predicted the oscillation of space and time known 
as gravitational waves (GW). Within a few years, with the operation of Advanced 
LIGO,^ Advanced VIRGO® and other sensitive interferometers, we anticipate the 
historic detection of GW. This will be a remarkable achievement and open a new 
window on the universe and on neutron rich matter. 

The first GW that are detected will likely come from the merger of two neutron 
stars. The rate of such mergers can be estimated from known binary systems.^ 
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During a merger the GW signal has a so called chirp form where the frequency rises 
as the two neutron stars spiral closer together. Deviations of this wave form from 
that expected for two point masses may allow one to deduce the equation of state of 
neutron rich matter and measure the radius of a neutron star r^s® Alternatively 
one may be able to observe the frequency of oscillations of the hyper-massive neutron 
star just before it collapses to a black hole. This frequency depends on the radius of 
the maximum mass neutron starP^ 1 

Continuous GW signals can also be detected, see for example.® Indeed Bild- 
stein and others® have speculated that some neutron stars in binary systems may 
radiate angular momentum in continuous GW at the same rate that angular momen- 
tum is gained from accretion. This would explain why the fastest observed neutron 
stars are only spinning at about half of the breakup rate. There are several very 
active ongoing and near future searches for continuous gravitational waves at LIGO, 
VIRGO and other detectors, see for example.® No signal has yet been detected. 
However, sensitive upper limits have been set. These limits constrain the shape of 
neutron stars. In some cases the star's elipticity e, which is that fractional difference 
in moments of inertia e = {I\ — I^)/^ is observed to be less than a part per million 
or even smaller. Here I±, I2, and I3 are the principle moments of inertia. 

An asymmetric mass on a rapidly rotating neutron star produces a time depen- 
dent mass quadrupole moment that radiates gravitational waves. However, one needs 
a way (strong stick) to hold the mass up. Magnetic fields can support mountains, see 
for example.® However, it may require large internal magnetic fields. Furthermore, 
if a star also has a large external dipole field, electromagnetic radiation may rapidly 
spin the star down and reduce the GW radiation. 

Alternatively, mountains can be supported by the solid neutron star crust. Re- 
cently we performed large scale MD simulations of the strength of neutron star 
crust.®'® A strong crust can support large deformations or "mountains" on neu- 
tron stars, see also,® that will radiate strong GW. How large can a neutron star 
mountain be before it collapses under the extreme gravity? This depends on the 
strength of the crust. We performed large scale MD simulations of crust breaking, 
where a sample was strained by moving top and bottom layers of frozen ions in op- 
posite directions.® These simulations involve up to 12 million ions and explore the 
effects of defects, impurities, and grain boundaries on the breaking stress. 

We find that neutron star crust is very strong because the high pressure prevents 
the formation of voids or fractures and because the long range coulomb interactions 
insure many redundant "bounds" between planes of ions. Neutron star crust is the 
strongest material known, according to our simulations. The breaking stress is 10 
billion times larger than that for steel. This is very promising for GW searches be- 
cause it shows that large mountains are possible, and these could produce detectable 
signals. 

To conclude this section, there is a great deal of interest in gravitational waves 
(GW) from neutron stars and there are many ongoing searches. One is interested in 
both burst sources, for example from neutron star mergers, and continuous sources 
from mountains or collective modes. Gravitational wave radiation depends on the 
equation of state of neutron rich matter. In addition, it can also depend on other 
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more detailed properties including the breaking strain of solid phases and the bulk 
and shear viscosities. 

§5. Conclusions: neutron rich matter 

Neutron rich matter is at the heart of many fundamental questions in Nuclear 
Physics and Astrophysics. What are the high density phases of QCD? Where did the 
chemical elements come from? What is the structure of many compact and energetic 
objects in the heavens, and what determines their electromagnetic, neutrino, and 
gravitational-wave radiations? Moreover, neutron rich matter is being studied with 
an extraordinary variety of new tools such as the Facility for Rare Isotope Beams 
(FRIB) and the Laser Interferometer Gravitational Wave Observatory (LIGO). 

We described the Lead Radius Experiment (PREX) that uses parity violating 
electron scattering to measure the neutron radius in 208 Pb. This has important 
implications for neutron stars and their crusts. We discussed X-ray observations of 
neutron star radii that also have important implications for neutron rich matter. 
Gravitational waves (GW) from sources such as neutron star mergers and rotating 
neutron star mountains open a new window on neutron rich matter. Using large scale 
molecular dynamics simulations, we found neutron star crust to be the strongest 
material known, some 10 billion times stronger than steel. It can support mountains 
on rotating neutron stars large enough to generate detectable gravitational waves. 

In conclusion, multi-messenger astronomy is based on the widely held belief that 
combining astronomical observations using photons, gravitational waves, and neutri- 
nos will fundamentally advance our knowledge of compact and energetic objects in 
the heavens. Compact objects such as neutron stars are, in fact, giant nuclei, even if 
they are an extraordinary 18 orders of magnitude larger than a 208 Pb nucleus. Nev- 
ertheless, both in the laboratory and in Astrophysics, these objects are made of the 
same neutrons, that undergo the same strong interactions, and have the same equa- 
tion of state. A measurement in one domain, be it Astrophysics or the laboratory, 
can have important implications in the other domain. Therefore we can generalize 
multi-messenger astronomy to multi-messenger observations of neutron rich matter. 
We believe that combing astronomical observations using photons, GW, and neu- 
trinos, with laboratory experiments on nuclei, heavy ion collisions, and radioactive 
beams will fundamentally advance our knowledge of the heavens, the dense phases 
of QCD, the origin of the elements, and of neutron rich matter. 
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